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Abstract 

The eigenvalues of the Corner Transfer Matrix Hamiltonian associated to the el- 
liptic R matrix of the eight vertex free fermion model are computed in the anisotropic 
case for magnetic field smaller than the critical value. An argument based on gen- 
erating functions is given, and the results are checked numerically. The spectrum 
consists of equally spaced levels. 



*PACS: 05.50.+q, 02.90.+p, 75.10Jm 







In recent times a great attention has been increasingly paid to the computation of 
correlation functions for integrable models in two dimensions. Following the approach of 
the Kyoto school it has been possible not only to diagonalize hamiltonians (as the XXZ 
model and its higher spin analogs || in the thermodynamic limit for the corresponding 
massive regimes) by means of the representation theory of quantum affine algebras, but 
also to extend the arguments to cases for which no quantum group like structure has 
been identified yet, namely the zero field eight-vertex (8V) model ||. A prominent role 
is played in this approach by the Corner Transfer Matrix Hamiltonian (CTMH), which 
in the quantum group invariant cases acts as a derivation in the quantum affine algebra 
(its physical interpretation is to be the lattice analog of the generator of Lorentz boosts 
JD). A salient feature of this operator is that its spectrum consists of equally spaced 
levels [[|. Remarkably there exists an 8V model (the elliptic free fermionic solution of 
the eight vertex Yang-Baxter equation || [7|, ||) for which a sort of affine quantum group 
structure CH q {2) has been identified, in which the R matrix acts as an intertwiner ||. 
The associated row-to-row hamiltonian is an (anisotropic) XY model in a magnetic field. 
This letter is a first step in the program of employing quantum group invariance to 
compute correlation functions in this 8V-like model. Namely it is shown here that the 
spectrum of the corresponding CTMH in the anisotropic regime 1 is given essentially by 
integer numbers, which makes it reasonable to try for it an interpretation as a derivation- 
like operator in CH q {2). This CTMH has already been studied in the context of the 
so called master symmetries [|H]] and its spectrum has been studied through the use of 
orthogonal polynomials (Tl|, [12|]. However in the latter references only the trigonometric 
(isotropic) and hyperbollic (disorder line) limits have been worked out, since the more 
general case leads to 5-term recurrence relations whose link to the 3-term recurrences 



characteristic of orthogonal polynomials [13fl is unclear. In the present letter the use of 



both a convenient parametrization of the model and an argument following Davies [13 1, 
dealing with generating functions enables us to identify the correct form of the eigenvalues 
of the CTMH also in a more general case. We compare the results to some numerical 
computations and then end with some comments and suggestions. 

The most general solution for the elliptic free fermionic 8V R matrix is, in the 
parametrization of ref. ||: 

R™(u) = l-e(«)e(^)e(^) R$(u) = e(u) - e(^i)e(^) 



1 for values of the magnetic field smaller than the critical one. See below. 



1 



R° l(u) = e(fo) - e(u)e(Vi) = e(V>i) - e(u)e(^ 2 ) 



(e(^ 1 )sn(^ 1 )) 1/2 (e(^)sn(^ 2 )) 1 / 2 (l - e(«))/sa(«/2) 
J2g(u) = -iA;(e(^i)sn(^ 1 )) 1/2 (e(^)sn(^)) 1 /2(i + e ( tt )) m (u/2) 

with e(u) the elliptic exponential: 

e(u) = cn(u) + i sn(u) (2) 

/c the elliptic modulus and cn(w), sa(u) the Jacobi elliptic functions of modulus k. In the 
sequel we will set 

i>\ = ?p2 = tp (3) 
since this is the interesting case from the hamiltonian point of view. Up to a normalization 
constant the row-to-row hamiltonian acting on the j-th and j + 1-th spins is 



. dR j j+1 (u) 



du 



;i + n^rrj., + (1 - V)ay j+1 + h{o* + o* i+x ) (4) 



n=0 



The cr's are the spin 1/2 Pauli matrices and 

r = k sn(^) , h = cn(-0) 



(5) 



For real ip and < k < 1, we stay in the Lieb-Schultz-Mattis (LSM) unit square [15], |T| 

o<r<i,o</i<i (6) 

We will take /i and k as our physical independent parameters, (see Figure 1) though some 
expressions will be left in terms of ip. A crucial point is the elliptic function dependence 
on this parameter shown in (p]). 



To define the CTM first construct a finite row of vertices by (see for instance [17]) 



Gf\u) = R n {u)R n „ l {u)---R j {u) 



(7) 



where Rj is the Boltzmann weight matrix correponding to the j-th vertex in the row. 
Then construct a wedge shaped region by piling up rows of vertices 

A N (u) = G[ N \u)Gi N \u) ■ --G^iu) (8) 



We shall consider the thermodynamic limit iV — > oo for which Baxter has shown [IS that 
the normalized (by its highest eigenvalue) CTM A^(u) is a one parameter commuting 
family written as 

A n (u) = exp(-uHcTAi) (9) 
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Expanding (§) for small u one finds for the CTMH 

oo 

F CT M = £j#;,;+i (10) 

This is the operator which we intend to diagonalize. Since the model is free fermionic we 
find it convenient to search for eigenvectors of (|TU|) which are linear combinations of the 



Jordan-Wigner fermions |19 



x.y 1 Hi V J 1 ( a f-i) x.y 



In terms of these 

oo 

Hqtm = ~2i£j {(1 + r)rfr/ +1 - (1 - + ^(r/r/ + t* +1 tJ +1 )} (12) 

i=i 

So look for states 

oo 
3=1 

such that 

[#ctm,V<0] = A,V(0 (14) 

Z is some index which labels eigenstates. We get 2 sets of equations for the numerical 
coefficients in (pi) (A, = -Aj/2, 7± = 1 ± V): 

7± (m - lJAj^.! - /i(2m - l)A± m + 7 T mA± m+1 = ±iA,AJ m (15) 

Note that if we substitute any of the (|T5| ) into the other we get a 5-term recurrence for 
either {A^}. Such relation readily becomes 3-term like if we set T 2 = 1 (Ising model 
|20|) or h = (doubled Ising model PD|). Also the trigonometric (k = 0) and the 
hyperbollic (k = 1, disorder line T 2 + h 2 = 1) limits of this model correspond in essence 
to three term recurrences which can be solved by means of orthogonal polynomials JTT . 
We want to solve for the anisotropic case (^), so define the functions 

oo 

A±(t) = ]T^± (16) 
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We get the differential equations 



( 7± t 2 - 2ht + 7 T 



dAf(t) 



dt 



+ (h- lT t- 1 )Af(t) = ±i~X l Ant) 



(17) 



Introducing the integrating factors 

f±(t) 

such that 



( 7± t 2 -2/it + 7 T )V2 
Af{t)=at{t)f ± {t) 



we are left with the system 



( 7+ t 2 - 2ht + 7 _) 1/2 ( 7 _t 2 - 2ht + 7+ ) 



1/2 



±\\iaf{t) 



Changing variables to 



dt 



( 7+ t 2 - 2ht + 7_)i/2( 7 _*2 - + 7+ )V2 



(19) 



(20) 



(21) 



the most general solution to (|2"0f ) becomes 



•l a e 



(0 - 
a, e 



X ') 



(22) 



where a® 2 are arbitrary integration constants. Due to the freedom in an overall constant 
we will be interested in the value of their ratio fax - So the whole issue resides 

in solving for the integral fl2i~D . Notice that in the generic case the integrand has four 
different poles (they appear as complex conjugate pairs) at 



cn(^) ± ik'sn(i[}) cn(-0) ± ik'sn(ip) 



1 + ksa(ip) 



1 — ksa(tp) 



(23) 



where k! is the conjugate elliptic modulus k 2 + k' 2 = 1. We recall the reader that these 
values coincide precisely with those of the central elements for quantum Clifford-Hopf 



CH q {2) invariant open chain Hamiltonians [pT |. 



An adequate change of variables is [22 



7+t 2 - 2ht + 7_ 
7_t 2 - 2ht + 7^ 
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(24) 



Notice that it breaks down in the trigonometric limit, which should be treated separately 

i±fe 

l=Ffe 

I fv+ dy 



in this approach. Being y± = the maximum and the minimum value of y respectively 
we get 



rv+ 
Jv 



2k'sn(ip) J y (y(y + -y)(y- y-)) 1/2 



(25) 



and then [23 



1 



tn 



-i 



\ 1/2 

1 + k \ t r_ — r j 



1 - k 



t r_i 



2k 1 / 2 
"' 1 + k 



(26) 



with r± = (1 ± sn(^)) 1 / 2 and tn _1 (w, k) an inverse Jacobian elliptic function of modulus 
k. Now we want that the functions flTB] ) generate normalizable states when N — > oo, 
so we require that they are free of singularities inside the unit circle \t\ = 1 This 
completely fixes the values of A; and with the result 



A/ 



2tt/ 
K'{k) 



|sn(V)| / = 0,1,2, 



(27) 



where K'(k) is the complete elliptic integral of the first kind of modulus k'. This spectrum 
reduces to its known values in the zero field {h — 0, ip — K(k)) [|T^] and disorder line [ P"T[ | 
cases. 

Finally we present some numerical confirmation of the above result for the spectrum. 
Since our CTMH is a quadratic form in fermionic operators we can try to diagonalize it 



by finding the eigenvalues Xf of the N x N LSM matrix fT5 

(A -B){A + B) 
where M is the following tridiagonal matrix 



MM 1 



(28) 



( h{0 + l) 7 _ 
7+ h(l + 2) 



M 



\ 










27+ 








27- 
/i(2 + 3) 
























(iV-2) 7+ h{N - 1 + JV - 2) (JV-1) 7 _ 
(JV-1) 7+ h(N-l) J 

(29) 



(Pq) is a pentadiagonal matrix as mentioned in and the eigenvalue equation for it 
leads to the same 5-term recurrence as flT5j). We have diagonalized it numerically for 
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k 


h 


C{ll), k)fhenr 

^-^ \T 5 / LlLCUf 


\-J \ lis ■ f\j 1 n n I'm 


lO" 2 


0.5 


0.908 


0.91 


0.1 


0.5 


1.47238 


1.47238 


1 


0.5 


3.46410 


3.46410 


0.5 


0.1 


2.89898 


2.89898 


0.5 


0.5 


2.52324 


2.52324 


0.5 


0.999 


0.1302 


0.1304 



Table 1: Numerical vs. theoretical values of C(ip,k). 



dimensions up to iV = 500, and found perfect agreement of the lowest eigenvalues with 
the formula ( p7|) for any values of k, h in the unit square @. In Fig. 2 we have plotted the 
lowest eigenvalues for fixed magnetic field h and let the anisotropy T vary with the elliptic 
modulus k; we find good fit of the theoretical coefficient C(ip, k) = 2'k\sh{'4j)\/ K'{k). This 
is also shown in Table 1. We have found poorer convergence in the cases of small tp or k 
(for fixed value of the other parameter), since they both correspond to approaching the 
isotropic limit T = (see Fig. 1) which is critical in this parametrization [pT|j . 

In conclusion we have found that the CTMH for the elliptic free fermionic eight vertex 
model has a spectrum given by integer values also in an anisotropic range of its param- 
eters. This leads us to think of its relation to some derivation-like operator for CH q {2) 
i0|| for its effective role as a derivation in the context of master symmetries), and 



sec 



makes it interesting to consider the problem of the computation of correlation functions 
in this 8V model using quantum group representation theory. This will be the subject 
of further study. It would be also interesting to complete the proof of orthogonality and 
completeness of the basis of eigenstates that has been generated here [0, which would 
involve the use of identities among elliptic functions, and to think of the pole structure of 
the quasimomentum-rapidity like transformation (|2lD in connection with quantum group 
invariance. Other interesting formal aspect is the relation (if any) of the 5-term recur- 
rence with the 3-term recurrences satisfied by orthogonal polynomials. The theory of 
(/-deformed orthogonal polynomials (see for instance ]24|]) might be useful in this. Finally 
since the trigonometric limit of the R matrix (|1]) is related to N = 2 SUSY (see 
for the more general CH q {D) case, D even) we are led to speculate whether the study 
of the algebraic properties of the CTMH we are considering could lead to the definition 
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of a lattice analog of an N = 2 superconformal algebra, much in the same way as [|l7j 
constructed a lattice analog of Virasoro (N = 0) for the free fermionic point of Baxter's 
solution of the zero field 8V model |18| (our h = case). 
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Figure Captions. 

Fig. 1. (h, k) parameter space. Our computation applies away from the two isotropic border 
lines h — 1, k — 0. 

Fig. 2. Lowest eigenvalues of MM 1 for fixed h = 0.5. Lines are guides for the eye. 
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